Radiation condition at infinity for the high-frequency 
Helmholtz equation: optimality of a non-refocusing 

criterion 

*Frangois Castella and ^Aurelien Klak 
March 2, 2013 

Abstract 

We consider the high frequency Helmholtz equation with a variable refraction index 
n 2 {x) (x G R d ), supplemented with a given high frequency source term supported near the 
origin x — 0. A small absorption parameter a £ > is added, which somehow prescribes 
a radiation condition at infinity for the considered Helmholtz equation. The semi-classical 
parameter is e > 0. We let e and a e go to zero simultaneaously. We study the question 
whether the indirectly prescribed radiation condition at infinity is satisfied uniformly along 
the asymptotic process e — > 0, or, in other words, whether the conveniently rescaled solution 
to the considered equation goes to the outgoing solution to the natural limiting Helmholtz 
equation. 

This question has been previously studied by the first autor in [ ]. In [ ], it is proved 
that the radiation condition is indeed satisfied uniformly in e, provided the refraction index 
satisfies a specific non-refocusing condition, a condition that is first pointed out in this 
reference. The non-refocusing condition requires, in essence, that the rays of geometric 
optics naturally associated with the high-frequency Helmholtz operator, and that are sent 
from the origin x — at time t — 0, should not refocus at some later time t > near the 
origin again. 

In the present text we show the optimality of the above mentionned non-refocusing 
condition, in the following sense. We exhibit a refraction index which does refocus the rays 
of geometric optics sent from the origin near the origin again, and, on the other hand, we 
completely compute the asymptotic behaviour of the solution to the associated Helmholtz 
equation: we show that the limiting solution does not satisfy the natural radiation condition 
at infinity. More precisely, we show that the limiting solution is a perturbation of the outgoing 
solution to the natural limiting Helmholtz equation, and that the perturbing term explicitly 
involves the contribution of the rays radiated from the origin which go back to the origin. 
This term is also conveniently modulated by a phase factor, which turns out to be the 
action along the above rays of the hamiltonian associated with the semiclassical Helmholtz 
equation. 
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1 Introduction 

1.1 General introduction 

In this article, we study the convergence as e approaches of w £ , solution to the following rescaled 
Helmholtz equation 

iea e w e {x) + ~w e {x) + n 2 (ex)w e (x) = S(x), x e R d (d>3). (1) 

Here a e is an absorption parameter, n 2 (x) is a space-dependent refraction index 1 and S(x) is a 
given and smooth source term. In the sequel, we assume the following: 

• The absorption parameter a £ satisfies 2 

a E > 0, a* — > 0. 

• The smooth refraction index n 2 (x) € C°°(M. d ) is a possibly long-range perturbation of a 
positive constant n 2 ^ > at infinity, namely, for some p > 0, we have 



VaeN d , 3C a , VxeR d 
where we denote as usual (x) := (1 + l^j 2 ) 1 / 2 . 



d«(n 2 (x)-nl) <C a {x)-"- a , (2) 



iflere and below we use the standard notation n 2 (x), a squared term, assuming in doing so that the corre- 
sponding term is everywhere non-negative. This is a harmless abuse of notation, since the refraction index n 2 (x) 
that is eventually chosen in our analysis is negative for certain values of x. The reader may safely skip this fact, 
since the Helmholtz equation also arises in the spectral analysis of Schrodinger operators, where the refraction 
index becomes E — V(x) where E is an energy and V(x) is a space-dependent potential, and the term E — V(x) 
may change sign in that context. 

2 The limiting case a e = + can be considered along our analysis, see below. 
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• The source term S(x) belongs to the Schwartz class 5(E). 

The question we raise is the following. Thanks to the absorption parameter a £ > in (1), 
the sequence of solutions w £ is uniquely defined (see below for the limiting case a £ = + ). On 
top of that, and as a consequence of specific homogeneous bounds obtained by Perthame and 
Vega in [14] (see [h] for extensions by Jecko and the first author, as well as [6]), it is clear that 
the sequence w £ is bounded in some weighted L 2 space, uniformly in e. Hence the sequence w £ 
possesses a limit (up to subsequences), say in the distribution sense, and the limit w = limw £ 
satisifies in the distribution sense the Helmholtz equation 

^w + n 2 (0)w = S, (3) 

where the variable coefficients refraction index n 2 (ex) in (1) has now coefficients frozen at the 
origin x = 0. 

Now, the difficulty is, the Helmholtz equation (3) does not have a uniquely defined solution. 
At least two distinct solutions exist, namely the outgoing solution, defined as 



w ut(x) : = lim ( iS + ~ + n 2 (0) ] S(x), 

5—>0+ \ Z J 



A, 



(4) 



and the incoming solution, defined similarly as Wi n — lim — i8 H \- n (0) S. Equiva- 

(5— >o+ V 2 J 

lently, the outgoing solution may be defined as the unique solution to the Helmholtz equation 
(3) which satisfies the so-called Sommerfeld radiation condition at infinity, namely 

A.V x iu out (x) +iV2n(Q)w out (x) = O ( t—^) , as \x\ — > +00. (5) 

\x\ VM / 

This formulation means that w out is required to oscillate like w ou t ~ exp (—iy/2n(0)\x\) /\x\ as 
I a; I — > 00. Similarly, the incoming solution satisfies the following radiation condition at infinity, 
namely (x/|x|) .V x Wi n — i"\/2n(0) iu, n = O (l/|x| 2 ), meaning that Wi n ~ exp (+i\^2 n(0)\x\) /\x\ 
as x —¥ 00. 

In that perspective, and due to the positive absorption parameter a £ > in (1), it is natural 
to expect that the previously defined sequence w £ goes to the outgoing solution w ou t to (3). 
This is the question we address here. 

It turns out that delicate analytical tools are needed to provide a clean understanding of 
the phenomena at hand, and to establish whether w £ ~ w out as e — > 0. The basic difficulty is a 
conflict between a local and a global phenomenon. On the one hand, the obvious fact that w £ goes 
to a solution to (3) is local: locally in x, i.e. in the distribution sense, the variable refraction index 
n 2 (ex) goes to the value n 2 (0) at the origin. On the other hand, the positive absorption parameter 
a £ > in (1) somehow asserts that w £ is an outgoing solution to A x w £ /2 + n 2 (ex) w £ — S, hence 
introducing the value at infinity = lim n(ex) = lim n(x), the solution w £ should roughly 

x— foo x— foo 

oscillate like w £ ~ exp (— i y/2 Uoo|x|) /\x\ at infinity. This is a global phenomenon. Now, all 
this is to be compared with the fact that w ou t oscillates like w ou t ~ exp (— i y2 n(0)|x|) /\x\ at 
infinity. Due to the fact that 7^ n(0), the radiation condition at infinity satisfied by w £ for any 
positive value e > is a priori incompatible with the radiation condition at infinity satisfied by 
the expected limit w ou t'. the radiation condition at infinity cannot be followed at once uniformly 



3 This assumption may be considerably relaxed at the price of some irrelevant technicalities. 
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in e, in any direct fashion (this is not in contradiction with the expected local convergence of w e 
towards w ou t-) 

Before going further, let us mention that the above question stems from a series of articles [1], 
[ ] about the high-frequency Hclmholtz equation (Equation (1) is a low-frequency equation) (see 
also [9] and [10] for similar considerations, in the case of a discontinuous refraction index, as well 
as [16] and [17] for the case of a variable absorption coefficient). These two papers investigate 
the high-frequency behaviour, in terms of semi-classical measures, of high-frequency Hclmholtz 
equations of the form 

iea E u E (x) + %- A x u e (x) + n 2 (x)u e (x) = J^S (^j (x 6 R d ). (6) 

The link between the low-frequency equation (1) that is the purpose of this article, and the 
high-frequency equation (6) is provided by the following basic observation: the function w e 
satisfies (1) if and only if the rescaled function 

U ^ X ) = ^dj2 W ^) ( 7 ) 

satisfies (6). In that picture, the main phenomenon to be described in (6) is the possibility of 
resonances between the high-frequency waves selected by the Hclmholtz operator e 2 A x /2+n(x), 
and the high-frequency waves carried by the rescaled source term e~ d l 2 S{x/e), both having the 
same wavelength e. Amongst others, it is established in [1], [3] that the semiclassical measure 
associated with u s can be completely computed provided w e indeed converges towards w ou t, this 
latter requirement being left as a conjecture in the cited papers. This is the motivation for the 
question we address here. 

In [4], the first positive convergence result w e — > w ou t is established. This results requires, 
amongst others, a specific and original non-refocusing condition on the refraction index n 2 {x) 
(called "transversality condition" in the original paper). This condition (see below for details) 
roughly asserts that the rays of geometric optics associated with the the semi-classical Hclmholtz 
operator e 2 A x /2 + n 2 (x) cannot focus at some positive time t > near the origin x = when 
issued from the origin at time t = 0. Later, X.-P. Wang and P. Zhang [19] proved a similar, 
positive result, using a so-called virial assumption which is stronger than the above non-refocusing 
condition. J.-F. Bony in [ ] establishes along quite different lines a positive result that is similar 
in spirit, requiring a weaker non-refocusing condition. 

The goal of the present text is to prove in some sense the optimality of the non-refocusing 
condition pointed out in [4]. 

We construct a refraction index n 2 (x) which violates the non-refocusing condition (rays of 
geometric optics issued from the origin do refocus close to the origin at some later time), and, by 
explicitly computing the asymptotic behaviour of w £ thanks to an appropriate amplitude/phase 
representation developped in [ 1 ] , we prove that 

w £ ~ w ou t + perturbation, 
e-s-0 v w / 

where the perturbation is computed as well. It explicitly involves the contribution of the rays 
issued from the origin which go back to the origin at some positive time, modulated by a phase 
factor that is the action, along these rays, of the hamiltonian associated with the high-frequency 
Helmholtz operator. 
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1.2 The non-refocusing condition 

As already mentionned, the asymptotic behaviour of w £ is dictated by that of the rescaled 
function u £ (x) = e~ d / 2 w £ {x/e). The function u £ is w E rescaled at the semi-classical scale, 
see (6) and (7). This is translated by the following identity, valid for any smooth test function 
4> € S(M. d ), namely 

where we denote as usual {w £ ,cf>) := / w £ (x) <fi* (x) dx, and * denotes complex conjugation. In 

JR d 

other words, the weak limit (w £ , (f>) of w £ can be computed as the weak limit at the semi-classical 
scale of u e , namely the limit of (u £ , e^ d ^ 2 (f>(x/e)) . This first observation is the main reason why 
semi-classical tools play a key role in our analysis. 

Besides, the asymptotic study of (1) is done here by transforming the problem into a time- 
dependent problem. This approach, introduced in [4], has been used since by J. F. -Bony ([2]) 
to study the Wigner measure associated to (6), or by J. Royer ([16]) when the absorption a £ 
depends oni. It consists in writing the solution w £ as the integral over the whole time of the 
propagator associated with iea £ + A x /2 + n 2 (ex), namely 

/■-f-OO 

w e (x)=i e~ a ' t e<^ +n2[£x) ) S{x)dt. (8) 



In the same way the outgoing solution can be written as 

r+oo 



w out (x): = i e lt (^+ n ' 2 ^) S(x)dt. 



In that picture, proving or disproving the convergence w £ ~ w ou t reduces to passing to the limit 
in the above time integral. 

Combining the two above observations, the basic first step of our analysis consists in writing, 
for any given test function </>, an in [4], 

(w £ , 4>) = (u e , e-^^x/e)) 

(U £ {t) S £ , <j> £ ) dt, (9) 

where we use the notation 

S £ (x) := -J^ S , and similarly <f> e (x) := <j) (^j , (10) 

where the semi-classical propagator associated with the semi-classical Hamiltonian e 2 A x /2 + 
n 2 (x) is 

t (e 2 



U e (t)=exp(i- ^Y A * + n ( x ))) ( n ) 

It is fairly clear on formula (9) that the asymptotics e — > in (w £ , 4>) is dominated on the one 
hand by the concentration of the rescaled test function cf> £ close to the origin at the semi-classical 
scale e, and on the other hand by the oscillations induced by the semi-classical propagator U £ (t) 
at the semi-classical scale s as well. The point is to measure the possible constructive intereference 
between both waves. 
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As standard in semiclassical analysis we define the semiclassical symbol 



h(x,0= l ^-n 2 (x), (12) 

associated with the semiclassical Schrodinger operator — ~ A x — n 2 (x). The semi-classical prop- 
agator U e (t) is known to roughly propagate the information along the rays of geometric optics, 
defined as the solutions to the Hamiltonian ODE associated with h, namely (see e.g. [8], [13], or 
[15]) 

-X{t, x, = S(t, x, 0, X(0, x, = x, 

< g (13) 

^ Q-Z(t,x,0 = V x n 2 {X{t,x,0), S(0,x,O=e 

It is clear as well that the integral J + °° ... in (9) carries most of its energy, semi-classically, over 
the zero energy level of h, defined as 

H Q := {(x,£) 6R M , s.t. h(x,£) = 0}. (14) 

In view of the integral (9) and of the above considerations, the following definitions are natural. 
The first definition is standard. 

Definition 1.1. [non-trapping condition] 

The refraction index n 2 is said non-trapping on the zero energy level whenever for each (x, £) € 
Hq, the associated trajectory (X(t,x,^),S(t,x,^)) satisfies 

lim \X(t,x,£)\ = +oo. 

t— 5- + 00 

When the refraction index is non-trapping, the rough idea is that any trajectory X(t, x, £) on 
the zero energy level leaves any given neighbourhood of the origin x = in finite time, making 
the above integral f Q + °° ... in (9) converge with respect to the bound t = +oo. 

The second definition comes from [ ] (this assumption is called "transversality condition" in 
the original text). 

Definition 1.2. [non-refocusing condition] 

We say that n 2 satisfies the non-refocusing condition if the refocusing set, defined as 

M := j(i,e,??) e]0,+ W [xl M s.t. ^ = n 2 (0), X(t,0,O = 0, S(t,0,O = v\ (15) 

is such that M is a submanifold of ]0, +oo[xM. 2d and M satisfies 

dimM < d - 1. 



When the non-refocusing condition is satisfied, the rough idea is that the trajectories X(t, 0, £) 
on the zero energy level issued from the origin x = at time t — cannot accumulate in any 
given neighbourhood of the origin x — at later times t > (this is encoded in the requirement 
on dimM). Technically speaking, an appropriate stationary phase argument in formula (9) 
allows to exploit in [4] the non-refocusing condition and to prove the weak convergence of w £ 
towards w out under this assumption. The main result in [4] is the following: when the refraction 
index is both non-trapping and satisfies the above non-refocusing condition, then w £ ~ w ou t as 
e — > weakly. 
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Recently, J.F. Bony in [2] shows the convergence of the Wigner measure associated with w e . 
He requires a geometrical assumption on the index of refraction that is in the similar spirit, yet 
weaker, than the above non-refocusing condition, namely 

meas n _! U e v / 2n 2 (0)S d - 1 ; 3t>0 X(t,0,£) = oj = 0, (16) 

where meas„_i is the Euclidian surface measure on v /2^2(0)S d " 1 and S^ 1 denotes the unit 
sphere in dimension d. Besides, inspired by [4], he constructs a refraction index which is both 
non-trapping and does not satisfy condition (16), and in that case he proves the non-uniqueness 
of the limiting of the Wigner measure. 

The goal of this paper is to construct a refraction index that is both non-trapping and violates 
the non-refocusing condition, and to establish in that case that w e goes weakly to a function of 
the form "wout+perturbation", for some explicitly computed and non-zero perturbation. To be 
more accurate, we construct below a refraction index for which the above refocusing manifold 

|„|2 

M = {(t,£_,r]) s.t. = n 2 (0), X(t,0,£) = 0, S(t, 0,£) = ?/} is smooth, yet has dimension 
dimAf = d — 1, a critical case, and we prove w e ~ "w ou t +perturbation " in that situation. 

1.3 Construction of the refraction index and statement of our main 
result 

Let us first examine the case of dimension d = 2. Let M s be a circular mirror centered at the 
origin. Any standard ray issued from the origin x = hits the mirror and goes back to the origin 
at some later time: refocusing occurs in a strong fashion. However all rays are trapped inside 
the circular mirror, leading to a trapping situation, in the sense of definition 1.1. To recover 
a non-trapping and refocusing situation, it is necessary to consider an angular aperture of the 
circular mirror, with total aperture < n. This is shown in figure 1: the circular mirror with 
total aperture < n provides a (non-smooth) non-trapping and refocusing refraction index. To 
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Figure 1: Spherical mirror in dimension 2 

transform the above paradigm into a smooth one, some regularizations need to be performed. 
The construction needs to be done in any dimension d > 2 as well. 
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Let us first introduce the hyperspherical coordinates (r, 9±, . . . , Qd-i) in dimension d > 2 

xi = rcos(0i), 

x 2 = r sin(^i) cos(# 2 ), 

a; 3 = r sin(0i) sin(# 2 ) cos(# 3 ), 

x d -i = rsm(0i)...sin(0 d _ 2 )cos(0d-i), 
a; d = rsin(0i)...sm(0d_ 2 )sm(0d-i), 

with 

#i G [0, 7r], 0j G [0, 27r] whenever j > 2 when g? > 3, and f?i G [— 7T, 7r] when d = 2. 
Next, we choose a fixed, smooth cut-off function x on IR such that 

X(*) = l, V|*|<1, x(*) = 0, V |*j > 2, x(t)>0, Vf G R. (17) 
We choose a radius R > and define the radial function 

f(x) = f(r) := X (2(r-R)), Vi = (r ) « ll ...,J i . 1 ). (18) 
We choose an angle (aperture) G [0, 7r/4[, and define the angular function 

g[x) = giOr) := X , Vx = Mb->M. (19) 

a smooth version of the angular aperture \8i\ < 8q. Finally, we choose two parameters > 
and A > such that 

nl < A. (20) 

We introduce the following 
Definition 1.3. [refraction index] 

We define the refraction index, retained in the whole subsequent analysis, as the following smooth 
version of the circular mirror with total aperture 9q < 7r/4 ; namely^ 

n 2 (x):^n 2 oc ~Xf(x)g(x) = n 2 oc ~Xf(r) 9 (9 1 ), \f x G R. (21) 

We are now in position to state our main result. Let (ei, . . . , e^) be the canonical basis of K d . 
Since the direction e\ is a symmetry axis for our refraction index, we introduce for later purposes 
the space Md(M) of square matrices of dimension d, we denote by 0<j(R) the space of orthogonal 
matrices, and we introduce the notation 

Q d ,i(K) := {A G O d (R), s.t. Aei = e x } . (22) 

The refraction index n 2 (x) in (21) is invariant under the action of Od,i(K). We last introduce a 
particular set of speeds, namely the set of initial speeds £ such that the zero energy trajectory 
X(t, 0, £) issued from the origin at time t = is reflected towards the origin at some later time 
t > 0. With the retained value of n 2 (x), we arrive at the definition 

4 The refraction index is negative in a bounded region of x. As already mentioned, we still use the abuse of 
notation consisting in using the squared of n. 
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Definition 1.4. [reflected rays] 

The reflection set Ig is defined as 

Ig = {Z:=(\Z\,6 1 ,...,e d _ 1 )€R d s.t. lG [-0 o ,+0o] and |f| = ^2^(0)} . 

Note that the (intuitive) fact that a velocity £ is such that X(t,0,£) hits the origin at some 
time t > if and only if £ £ Ig , is proved later (see section 2.2). 

Our main result in this text is the 
Theorem 1.5. [Main Result] 

Let n 2 be the refraction index defined in (21). Assume the aperture 9q < 7r/4 and the radius 
R > satisfy the smallness condition 

l-cos(20 o ) < (23) 
Assume d > 3. Then, the following holds: 

i) The index n 2 is non-trapping on the zero-energy level Hq — {(x, £) s.t. |£| 2 /2 — n 2 (x) = 0}. 

ii) The refocusmg set M = {(t,£,?j) s.t. |?/| 2 = 2n 2 (0), -X"(f,0,£) = 0, E(i, 0, f ) = ??} fsee fi5jj 
is a smooth submanifold o/]0, +oo[xK 2d , boundary, and its dimension has the critical value 

dim(Af) = d- 1. 

in) Assume the source term S satisfies S G 5(R d ). Then, we have 

V(f>eS{E. d ), (w e - (w out + L e ) , 4>) — > 0, 
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where the distribution L e is defined for any <j> G 5(K d ) through 

(Le,<l>) = 

r,,,,, / ..-x,»j •-- I I ^ E i s iMt +n 2 (x( s ,o,0) ) ds) S(0^(-Z)da ea (0. (24) 




Here dae denote the natural Euclidean surface measure on Ig (see definition 1-4), the return 
time Tr > is the unique time^ such that for any £ € Ig we have X(Tr,0,£) = 0, and the 
constant C n 2 d ^ can be explicitly computed and depends only on the index n 2 and on the 
dimension d. 

Remark. The condition (23) is technical, and requires the aperture 9q to be small: it ensures 
the trajectories cannot be trapped by the refraction index. 

Remark. Note in passing that the constraint d > 3, which is also needed in reference [4], comes 
from a stationary phase argument. This constraint on the dimension is standard in the analysis 
of Schrddinger-like operators. It comes from the fact that the dispersion induced by the free 
Schrddinger operator acts like t~ d l 2 , a factor that is integrable close to t = +oo whenever d > 3. 

Remark. Let £o := {y2 n(0), 0, . . . , 0). The distribution L e can as well be written as 

(Le , 0) = 

C„V ex P I ~ / I '" (a '°' go)|2 +n 2 (X( S ,(Uo)) ) ds ) I / S(0**(-0d*flo(O 

This formulation illustrates in a clearer way the fact that if the source S radiates towards the 
mirror, then w £ converges towards a non-trivial perturbation ofw out . 

Note in passing that in the present counter-example, as in the paper by J.-F. Bony [ j , only 
subsequences of w £ converge, due to the above oscillatory factor exp(i const, je). 

Remark. In the chosen hyperspherical coordinates, the Euclidean measure da$ (£) coincides 
with do~g (l;) = n(0) d ~ 1 da(9\, . . . , Od-i), where da(9i, . . . , 9d-i) denotes the standard euclidean 
surface measure on the unit sphere S d_1 . 





1.4 Preliminary reduction of the proof 

Our main result contains three distinct statements. Items (i) and (ii) are of geometric nature, and 
merely concern the behavious of the classical trajectories associated with the retained refraction 
index. Their proof is performed in sections 2.1 and 2.2, respectively. Item (iii) is the main item, 
and concerns the asymptotic analysis of w e . Since our analysis heavily relies on tools previously 
developped in [ ]. we briefly recall here some of these tools and indicate how the analysis of w e 
can be reduced to a simpler sub-problem. We postpone the analysis of the reduced subproblcm, 
hence of item (iii) of our main result, to section 3 below. 

As already indicated, given a smooth test function <f>, we start from the formulation 

(w s ,(f>) = - e- a ^(U E (t)S e ,cl> e )dt 
£ Jo 

5 The fact that all these quantities exist and are well defined is part of the Theorem, and is proved in section 
2.2. 
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(See above for the notation) . The next step consists in splitting the above time integral into four 
time scales, namely very small, small, moderate, and large time scales. To do so, we take one 
small parameter 6 > and two large parameters Tq > and T\ > 0, and split the above time 
integral into the four zones 

< t < T e, T e<t<9, 9<t<T u T x < t < +oo (6 < 1, T ,T r > 1). 

Technically, we use a smooth splitting, based on the already used cut-off function x ( see (17)). 
Besides, we also distinguish between the contribution of zero and non-zero energies, namely 
taking a small parameter 5 > 0, we write, in the sense of functional caculus for self-adjoint 
operators, the identity 

2 

1 = X s(H e ) + (1 - xs){H e ), where H e := E — A x + n 2 (x), and X s(s) := X (|) (fl G K, <5 « 1). 

The main intermediate result of the present subsection is the following 

Proposition 1.6. [Main intermediate result] 

Take a test function cf> 6 l S(M' i ). Define uT e as 

(w e ,cj>) := i ^ (1 - X ) e-^ <l7 e (t) x«(F e ) ft, 0s) dt. 

Then, there is a large T± > sitc/i that for any small 8 > 0, and any small 9 > 0, f/iere exists a 
constant Ce,s > suc/i £/ia£ /or any small e > 0, we /law 

(to e - (w «t + 53^)) , 4> 

This result roughly asserts that w e is asymptotic to w out + u% as e — > 0, up to carefully 
choosing the various parameters Tq, T\, etc. Hence the proof of item (iii) of our main result 
essentially reduces to proving that uT e ~ L £ as e — > 0. 

Proof of Proposition 1.6. 

The proof is obtained by gathering the statements of Proposition 1.7, Proposition 1.8, Propo- 
sition 1.9, Proposition 1.10 below. ■ 

The remainder part of this paragraph is devoted to a brief idea of proof of the above auxiliary 
Propositions that lead to Proposition 1.6. 



<Cg, 5 



1 



d/2-1 



t; 
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• Contribution of very small times < t < Tq e. 

r+oo 

The contribution of very small times to (w £ , cf>) = ie^ 1 I e~ aet (U e (t)S £ , <f> £ ) dt, is 

Jo 

i f 2Toe / t \ 

It is the main contribution to w £ , provided To is large enough. Indeed, we have the following 
fact, whose proof is based on a simple weak convergence argument. 
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Proposition 1.7. (See [ ]). Let n 2 (x) be any bounded and continuous refraction index. Then, 
if S and (j> belong to S(Mr), we have 



(i) For all time T > 0, 



2T e 



T e 



e- a ^{U e (t)Se,^)dt^i X 



2T„ 



dt. 



(ii) There exists C c i > which only depends on the dimension such that 



X f — ) < exp(it(A x /2 + n 2 (Q)))S, (j>)dt)- (w ouU , 



< 



C d 



d/2-l ' 



T, 



• Contribution of small, up to large times, away from the zero-energy level. 

r+oo 

The contribution to (w e , <f) = i e^ 1 / e _Qe< (U e (t)S £ , </> e ) dt that is associated with small, 

Jo 

up to large times, away from the zero-energy level, is 



s (i-x) 



Toe 



TnS 



((l-Xs)(.H e )U e (t)S e ,<f> c )dt. 



It is seen to be small, using a non-stationary phase argument in time, see [4] (this is the reason 
for the previous cut-off close to the initial time t — 0, where integrations by parts in time are 
forbidden). Indeed, we have the 

Proposition 1.8. (See [ ]). Let n 2 be any long-range refraction index. Let S and <f> belong 
to L 2 (R d ). Then there exists a constant Cs > 0, which only depends on S > 0, such that for any 
small e > and any Tq > 0, we have 



(i-x). T . 

Toe \ 1 0£ 



((l- X s(H e ))U E (t)S E ,<f> £ (t))dt 



• Contribution of large times, near the zero-energy level. 



The contribution to (w E ,<j>) =ie 
times, close to the zero-energy level, is 

f + OO 



r+oo 

1 / e~ act (U e (t)S e , <f> e ) dt that is associated with large 
Jo 



(Xs(H e )U s (t)S £ ,<p e )dt. 



It is seen to be of order 0(e N ), for all N € N, see [4]. Indeed, the semiclassical support of 
Xs(H E )U E (t)S e goes to infinity in the x direction at speed of the order 1 (i.e. the semi-classical 
support lies in a region that is at distance of order t from the origin - this uses an argument due to 
Wang, see [18]), while the semi-classical support of <fi e remains close to the origin. This argument 
relies on the fact that for T\ large enough, the semiclassical supports of the two functions are 
disconnected, which in turn uses the non-trapping behaviour of the refraction index. We arrive 
at 
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Proposition 1.9. (See [ ]). Let n 2 be any long-range refraction index that is non-trapping. 
Let S and be in S(M. d ). Then there exist Sq > and Ti(Sq) > such that for all time T\ > Ti(Sq) 
and any < 5 < 5q, there exists a constant Cs such that 



-foo 



< C s e. 



• Contribution of small times near the zero-energy level 

f+OO 



/•-l-oo 

The contribution to (w £ ,<ft) = ie^ 1 / e~ aEt (U £ (t)S £ ,(j) £ ) dt that is associated with small 

Jo 

times, close to the zero-energy level, is 

e"^ (1 ~ x) (^) X (U E (t)Xs{H e )S e , cf> e )dt. 



Unlike in the previous case, the semiclassical supports of U £ (t)xs(H £ )S £ and <p £ may intersect 
for these values of time t. The whole point in [4] lies, roughly speaking, in proving a dispersion 
estimate. The key is to prove that the variable coefficients Schrodinger propagator U £ (t) has 
the same dispersive properties than the free Schrodinger propagator, corresponding to the case 
when n 2 = 0, at least for small values of t such that < t < 9 (for later times, the semiclassical 
support of U £ (t) S £ is close to the classical trajectories (X(t), S(i)), trajectories which in turn may 
come back close to the origin and contradict any dispersion effect). Indeed, for small times, the 
trajectory (X(t), S(i)) is close to its first order expansion in time, which is the key to obtaining 
dispersive effects similar to the one at hand in the free case. Technically speaking, the proof relies 
on establishing that the propagator U £ (t) behaves like the free Schrodinger propagator for small 
times, a propagator whose symbol is exp(ii|£| 2 / e) , and which in turn has size (e/t) d/2 thanks to 
a stationary phase argument. 

To obtain the desired statement, a wave packet approach is actually introduced, which 
strongly uses the work by Combescure and Robert ([ ]). It allows to compute explicitly the 
propagator U £ (t) S £ , using the Hamiltonian flow and related, linearized, quantities, to obtain a 
representation of the form 



e 

T e 



(1 - X) (^A X (j) (U £ (i) XS {H £ )S £ , £ ) dt 

e*^ a ' a N (t,X)dtdX + Og :S (e N ) , (25) 



1 



r(5d+2)/2 



T e JM. 6 



where X = (q,p,x, 2/, ^, 17) G K 6d , where TV is a possibly large integer, and the remainder term 
Og : s (e N ) is upper bounded by Cgj £ N for some Ce,s > independent of e, which depends on 
the chosen 9 > and S > 0. Note that the amplitude ajv is defined in (38) below, while the 
complex phase function ip is defined in (37) below. We refer to section 3 for details about the 
representation formula (25), which is a key ingredient in our proof of the main theorem. 

With this representation at hand, we arrive at the 

Proposition 1.10. (See [ ]). Let n 2 be any long-range potential which is non-trapping. For 
9 and 5 small enough, there exists Cg > and Cg^s > such that for all s < 1 we have 

1 (l-x(^))^ t (t/e(^( J ffeR,0s)^<^ T + C M£ . (26) 
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2 Properties of the refraction index 
2.1 Non-trapping behaviour 

The goal of this subsection is to prove item (i) of our main Theorem 1.5. 

We prove that the chosen refraction index n 2 (x) = n 2 ^ — Xf(r)g(6x) in (21) is non-trapping 
on the zero-energy level H = {(x,£) € R 2d , s.t. £ 2 /2 = n 2 (x)}. 

We first observe that the zero energy level has the more explicit value 

H = eR 2d , a.t.x = (r,6i,...,0 d -i), y = - A/(r) 5 (^)| . 

We readily define the following two regions. The first one is usually called the classically forbidden 
region: any trajectory living on the zero-energy level cannot reach the set J5g. The second one is 
sometimes called here the bump of the refraction index: it is the region where the refraction index 
actually varies with x. Outside this region, the refraction index is constant and the Hamiltonian 
trajectories associated with h(x,^) = |£| 2 /2 + n 2 {x) are straight lines. 

Definition 2.1. (i) We denote by Bq the set (classically forbidden region) 

B % := {x £ R d , s.t. n 2 (x) < 0} = {x = (r, X , . . . , 9^), s.t. n 2 ^ < A/(r) 5 (^)} . 

(ii) We denote by B p the set (bump) 

B p :={x=(r,6 1 ,,...,0 d - 1 ), s.t. R-l<r<R+l, \6 X \ < 26 } . 

Remark. From the definition of B§ and the two functions f(r) — x(2(?" — R)) an d ff(^i) = 
x(Qi/&o) it * s clear that there exists /i (E]l,2[ such that 

B $ C {R- I <r<R+^, </i^}. (27) 

It suffices to take [i such that 

Th 

0< X (m) < 77= (hence fie]l,2[). (28) 

Our main step lies in proving the following escape estimate 

Lemma 2.2. Select the refraction index n 2 (x) as in (21) and assume condition (23) is fulfilled, 
namely 1 — cos(2#o) < 1 / (2i?) . Take a Hamiltonian trajectory X(t, x,t;) = X(t) living on the 
zero-energy level and define Xg := (R, 0, . . . , 0) in Cartesian coordinates. 
Then, there exists a > 0, as well as /3 € R and 7 € M, such that 

Vt>0, \X(t)-x \ 2 > at 2 + fit + 7. 

An immediate corollary of the above Lemma is 

Corollary 2.3. Assume condition (23) is fulfilled, namely 1 — cos(2#o) < l/(2i?)- Then the 
refraction index n 2 (x) in (21) is non-trapping on the zero-energy level. 

Proof of Corollary 2.3. 

Apply the preceding lemma and let t — ¥ +00. M 

Proof of Lemma 2.2. 
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MB, 



R-l/2 R R+l/2 



R + l el 



Figure 3: Bump of refraction index, and classically forbiden region 



First step. We compute the second derivative of \X(t) — x \ and get 



dX 
~dt 



(*) 



= (Vn 2 (X(t)),X(t)-x ) + 



dX 
~dt 



(*) 



(Vn 2 (X(t)),X(t)-x )+n 2 (X(t)), 



(29) 



where we have used the fact that the Hamiltonian trajectory (X(t), £(t)) belongs to H . Letting 
X(t) = ru r in hyperspherical coordinates and x — (R, 0, ...,0) in Cartesian coordinates, we 
obtain on the other hand 

(Wn 2 (X(t)),X(t) - x ) = l—\f'{r)g{8{)u T - X^g'^ue^ritr - RSi) , 



= F r (r,0 1 ) + F e (r,0 1 ), 



where 



F r (r, 9 1 ) = -A/'(r)2(0i) (r - ficos(^)) , F e (r, 9 1 ) = —\—f(r)g'{6\) sin(0i). (30) 

r 



Eventually we have 



^ \X{t) -x Q \ 2 = F r (r, 9,) + F e (r, 9,) + n 2 {X{t)). 



(31) 



Therefore, the lemma is proved once we establish the existence of a > such that 

F r (x) + Fg{x) + n 2 (x) >a>0 
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whenever x G H x Hq = R d \ (where H x denotes the projection (a;, x from R 2d to Mr). 
We readily notice that n 2 and Fg are clearly non-negative function on the whole of W 1 . 

• Step two: non-negativity of F r . First, on M. d \ B p , the function F r is zero, hence non- 
negative. In the same way on B p n {R — 1/2 < r < i? + 1/2}, we have /' = 0, hence F r = > 0. 
There remains to study the non-negativity of F r on the two sets {R — 1 < r < R — 1/2, |6>i| < 28q} 
and{i?+l/2<r <R+l, \0i\ < 29 }. 

On {R - 1 < r < i? - 1/2, |0 X | < 26> }, we have 

r - i?cos(6»i) < J? - ~ - i?cos(26» ) = i?(l - cos(2(9 )) - - < 0, 

thanks to our assumption (23). Since /' > on {R — 1 < r < R — 1/2}, we get F r > 
on {R — 1 < r < R— 1/2, |#i| < 2#o}- A similar computation proves that F r > on the set 
{R + l/2<r<R + l, \6 1 \< 26 }. 

We have obtained that F r > on the whole of Mr. 

• Step three: decomposition of R d . We have just proved that F r (x) + Fg(x) + n 2 (x) > 
for all x € R d . We now wish to obtain a positive lower bound for x ^ B%. The argument 
relies on the fact that the refraction index n 2 is positive away from the boundary where 
dB$ :— {(r, #i, . . . , #d-i)> f(r)g(0i) — n^/A}, while the term F r + Fg stemming from the 
gradient of the refraction index in (31) is positive close to the boundary dB$. This is the reason 
for the decomposition we now introduce. 

We define the set (piece of ring) 

C a ,0 :={R-a<r<R + a 7 -0 < Q x < (3} . 
We know from the remark after Definition 2.1 that there exist /i s]l, 2[ such that 

Bq C C R+ n/2 jl/t g . 

We therefore decompose 

Rd \ B <& = ( R \ Ci?+M/2,M0o) U { C R+h/2,^8 \ -B0) ■ 

We readily observe that, by construction of \x (namely x(m) 2 G]0; n L/-M ~ see (28)), for any 
iGI^ Cfl+^/2,^e 1 we nave the lower bound 

n 2 (x) = - Xf(r)g(9 x ) > - A X (,i) 2 =: c n , > 0, 

There only remains to prove the existence of cy > such that F r + Fg > cv on C R+tl /2,fj,6 \ ^0- 

• Step four: positive lower bound for F r + Fg on C R+fl /2.n0 a \ B<b- Take v €]1,2[ such 
that 

7=x < x{v) < 1 

where \ is the truncation function defined in (17). With this choice of v, we clearly have, 
whenever x € C R+u / 2 , v g , the relation n 2 (x) = n 2 ^ - \x(2{r - R))x(6x/0 ) < n 2 ^ - A%(^) 2 < 0, 
hence 



1G 



Therefore, it is enough to obtain a lower bound on F r + Fg on the set Cn +fl /2,^8 \ Cr.+v/2,v8< 
To this end, we decompose (see Figure 4) 

CR+H/2,(t6 \ Cr+v/2,v6 CZ^U^UZjUZj, with 

Z\ := {R - n/2 < r < R - v/2, \Q X \ < vO Q } , 

Z 2 r := {R + v/2 < r < R + fi/2, <!/$)}, 

Z] := {R - n/2 <r<R + fi/2, < 6 l < -v6 } , 

Z 2 e := {R — fi/2 < r < R + /j,/2, v9 < 6 X < ii6 } ■ 




Figure 4: Zone of study 



On Z\. We use the structural hypothesis (23) to get 

F r {z) = -\f{r)g{6 x ){r - i?cos(^)) > -A/'(r) 5 (^)(i? - V - - Rcos(26 )) 

> A/'(r) 5 (^i)^ > X(u - 1) ( mm x'(s)) (mm X {s)) =: ci > 0. (32) 

A similar proof establishes that, whenever x £ Z 2 we have 

F r (x) > \{v - 1) ( min \-x'{ s )] I ( mm x( s ) ) =■ c 2 > 0. 
\se[v,n] J \\s\<is J 

On Zg. The important term is now Fg. We have 

F e (x) = -X^f(r)g'(e 1 )am(e 1 ) > Xjf(r)g'(e 1 )sm(ue ) 

* X e Q (R + ,/2) fe x(s) ) C 6[ T-,] X ' (s) ) =: C3 > a 
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A similar argument establishes that, whenever i 6 Z e wc have 

Fe(x) > A t— — f min J ( min [-rf (s)] ) =: c 4 > 0. 

Gathering all estimates, there exists a positive constant cv > such that 
Vx e C R+li / 2 ,po \ C R+I//2 , u e , F r (x) + F g (x) > c v > 0. 

• Step five: end of the proof. Putting all estimates together, we obtain 

Vx G ll x Ho = R d \ B $ , F r {x) + F g (x) +n 2 (x) > min(c n 2,c v ) =: a > 0. 
The lemma is proved. ■ 

2.2 Refocusing Set 

The goal of this subsection is to establish part (ii) of our main Theorem 1.5. 
Our main result is 

Proposition 2.4. Letn 2 be the potential defined in (21). Assume the structural hypothesis (23) 
is fulfilled, namely 1 — cos(2# ) < 1/(2R). Then, the refocusing set defined in Definition 1.2 as 

M=^(t^,rj)e}0,+oo[xR 2d s.t, = n 2 (0), X(t, 0, Q = 0, E(t, 0, = V 

satisfies 

M = {(T H , S .(.(=-f|=Mi,..,U r= ^2n 2 (0), |^| < 9 } , 

where Tr > is the unique positive time such that X(Tr, 0, (•y^2n 2 (0), . . . , 0)) = 0. 

Proof of Proposition 2.4. 

Consider a trajectory X(t, 0, £) = X(t) on the zero energy level, with ^ = (r, 6* 1; . . . , fl^-i) m 
hyperspherical coordinates. 

If 1 6*i | > 26*o, it is clear that X(t) is a straight line which never enters B p , and the equation 
X(t, 0, £) = with t > has no solution. 

We need to understand the geometry when the trajectory reaches B p , i.e. when \9%\ < 29q. 

We prove below that two cases occur. If \9\\ < 9q, the trajectory remains along a line, and 
it is reflected by the refraction index towards the origin. If 9q < |#i| < 29q, the force acting on 
the trajectory has a non-vanishing component in the orthoradial direction, which prevents the 
trajectory to go back to the origin. The proposition follows. 

Let us come to a proof. 

• First case: \9i \ < 9q. 

Consider the trajectory Y{t) defined in hyperspherical coordinates as 

Y(t) = (r(t),9 u ...,9 d ^), 
with r(t) solution to the ordinary equation r" = —Xf'(r) with initial data 

r(0) = 0, r'(0) = ^2^(0). 
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Then, (Y(t),Y'(t)) satisfies the Hamiltonian ODE (13) associated with h(x,£) = |£| 2 /2 + n 2 (x). 
Since Y(0) = Z(0) = 0, and Y'(0) = X'(0) = f, uniqueness provides = F(t) for all i. The 
trajectory is radial. 

It is clear that the radial trajectory t H> r(t) reaches the region {R — l<r<i?+l}at time 
t e = (R- 1)/|£| = (R- l)/y / 2n 5 (0) > 0, where £ e = inf {i > 0, X(t) G B p } . Now, according to 
Corollary 2.3, the trajectory r(t) necessarily leaves the region {R — 1 < r < R + 1} at some later 
time t s > t e , where t s — inf {t > t e , X(t) £ B p } . The trajectory can either leave the bump at 
r = R—l or at r = R+l. The case r = R+l is forbidden, for in the contrary case, using continuity, 
there would exist a time t c such that r(t c ) — R, hence X(t c ) € B$, which is not allowed. 
Therefore, the trajectory leaves the bump B p at X(t s ) where |X(i s )| = r(t s ) = R — 1. Energy 
conservation, together with the fact that the trajectory is radial, implies that X'(t s ) = — £. 
Therefore, the trajectory for later times t > t s is a straight line with constant speed — £. We 
deduce that there exists a unique Tr > t s such that X(Tjj,0,^) = 0, and we have as desired 
3(2^,0,0 = -6 

• Second case: do < \6i\ < 28o. 

We first assume that d = 2, and next generalize the argument to d > 3 using the symmetries 
of the system. To fix the ideas, we assume in the following that 0q < 9\ < 29o, the proof being 
the same when 9\ has the opposite sign. 
* In dimension d = 2. 

Let t e = (R— 1)/|£| be the time when the trajectory enters B p , as in the preceding case. 

On the one hand, since the velocity 2(i e ) is radial and satisfies S(i e ) = |£| u r , there is an e > 
such that R — 1 < < i? + 1 whenever t E]t e ,t e + e]. On the other hand, by assumption we 

have 9\{t e ) = 9\ e]9o,29o[, and continuity implies there is an e > such that 9q < 9\{t) < 29q 
whenever t € [t e , t e + e]. Hence we may define 

t s := sup{t > t e , s.t. W € [t e ,t], 9 X (t') €] 9 , 29 [ and X(t') ^ 0.}. 

Now, Hamilton's equations of motion (13) can be written in polar coordinates as 

r»-r{0' 1 ) 2 = -Xf'(r)g(6 1 ), 
2r>9' 1 +r9'{ = -\£&g'(9 1 ). 

Examining the second equation, we have (r 2 9' 1 )' = 2rr'9[ + r 2 9'{ — —Xf(r)g(6i), and we get 
whenever r(t) =/= 0, 

e[(t) = -~~J* f(r(s))g'(9 1 (s))ds. (33) 

Therefore, since f(r) > for any r > while f(r) > whenever R — 1 < r < R + 1, and since 
g'{9i) < when 9 < 9\ < 26 , while g'(9 1 ) < when 9 a < 9\ < 29 we get, with the above 
definitions and observations, 

9[(t)>0 Vt€]te,t s }. 

With this observation at hand, two cases may occur. 

If t s = +oo, there is nothing to prove, for by definition of t s , we have X(t) ^0 whenever 
< t < t s — +oo. 

In the case t s < +oo, we already know X(t) ^ whenever < t < t s . Besides, since 
9[(t) > whenever < t < t s , it is clear that the case X(t s ) = is impossible (for in that case 
the trajectory would be a straight line passing through the origin on some interval [t*,t s ], in 
contradiction with 9[(t) > on [t*,f s ]), hence 9x(t s ) = 29q and 9' 1 (t s ) > 0. For that reason, the 
trajectory X(t) for times t > t s is a straight line with constant velocity, which lies entirely in 
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the set 26>o < Oi < 28q + n. In particular, since 0'i(t s ) > 0, the trajectory cannot be radial and 
we have X(t) 7^ whenever t > t s in that case. This concludes the proof. 

* In dimension d > 3. 

We use the invariance of n 2 under the action of 0,j,i(K). 

Take £, E R d such that |£| = y/2n 2 {0). Write £ = (■ y /2n 2 (0),8 1 , . . . , in hyperspherical 

coordinates. There exists a matrix € ©d.i(R) such that £ = (-\/2n 2 (0), 9i, 0, ... ,0). On 
the other hand, denote by (r(i), 6*i(i)) the solution of Hamilton's equations of motion (13) with 
initial data (■ v /2n 2 (0), 6>i) in dimension 2. We set Y(t) = A^(r(t), 9 x {t), . . . , 0). Then F(t) 
satisfies Hamilton's equations of motion (13), with initial data Y(Q) — 0, y'(0) = £. Uniqueness 
provides Y(t) = X(t) for any t > 0. This, combined with the previous step, provides X(t) ^ 
for any t > 0. ■ 

3 Convergence proof 

The goal of this section is to prove item (iii) of our main Theorem 1.5. 

The proof is performed in a number of steps. We begin by defining some necessary notation. 

3.1 The linearized hamiltonian flow 

Let ip(t,x,£) = (X(t, x, £),S(i, x, £)) denote the flow associated with Hamilton's equations of 
motion (13). The linearized flow, written F(t, x, £) below, is 




where A(t), B(t), C{t), D(t) are by definition 



A(t,x,0 



C(t,x,0 



Dx 
Dx 



D(t,x,0 



DC 
DZ(t,x,Q 



The linearisation of (13) leads to 




Id, 



(34) 



as well as 




Id, 



0. 



(35) 



Finally, we define for later purposes the matrix T(t, x,£) as 



r(t, x, = (C(t, x,0 + i D(t, x, 0) • (A(t, x, + iB(t, x, £)) 



-l 



(36) 
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3.2 A wave packet approach: preparing for a stationary phase argu- 
ment 

The intermediate result in Proposition 1.6 establishes roughly that (w £ ,(f>) ~ (w ou t + w £ ,(f>} as 
e — > 0. Therefore, item (iii) of our main Theorem reduces to proving (w^, 4>) ~ (L e , 4>) as s — > 0. 
Therefore, this preliminary paragraph is devoted to express the quantity 



1 



(.<?:■<>) = - £ r ( i - \ (;))<• "•' :r : (/) X ,r/./ ; !,s' : .o : },//. 



as an appropriate oscillatory integral. Our approach uses the technique developped in [4], which 
in turn strongly uses a wave packet theorem due to M. Combescure and D. Robert (see [7]). We 
skip here the details of the proof, referring to [4]. 
The main result in this paragraph is the following 



Proposition 3.1. (See [ ]) Whenever X = (q,p, x, £,y,r)) € R andt £ K, define the complex 
phase 

^(t, X) := J + n 2 (q s )^j ds - p.(x -q)+ Pt .(y- q t ) 

+ x4 - y.v + i ^ x ^ q)2 + ir t (y - qt ).(y - qt), (37) 

where q t :— X(t,q,p), pt ■= S(t, q,p), and T t := T(t,q,p). Select an integer N G N. Select two 
truncation functions Xo(<7>p) an d Xi{ x >y) both lying in C^°(M. 2d ), and such that 

supp X o(q,p) C < 26} U {\\p\ 2 /2 - n 2 {q)\ < 28} , 
Xo(q,p) = l on {\q\<36/2}u{\\p\ 2 /2-n 2 (q)\<35/2}, 
Xi(x,y) = 1 close to (0,0). 

Define the amplitude 

a N (t,X) -e-^il-x) (I) Sm*(v)xv(q,p)xi(x,y)PN [t,q,p,^^j , (38) 
where P]\r(t,q,p, z) satisfies 

P N (t,q,p,x) := ^ JI det(A{t,q,p) + iB{t,q,p))- 1 / 2 Q N (t,q,p,x), (39) 

and the square root det(A(t, q,p) + iB(t,q,p)) c is defined by continuously following the argu- 
ment of the relevant complex number, starting from the value det(A(0,q,p) + iB(0,q,p) = 1 at 
time t — 0, while Q^(t,q,p,x) is a polynomial in the variable x G K d ; whose coefficients vary 
smoothly with (t,q,p), and e, and which satisfies 

Q N (t,q,p,x) =l + 0(Ve) 

in the relevant topology. More precisely, we have 

Q N (t,q,p,x) = 1+ ^2 £ ^ J Pk,j(t,q,P,x), 

(k,j)ei N (40) 
In = {1 < 3 < - 1, 1 < k - 2j < 2N - 1, k > 3j} , 
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where each pk.j has at most degree k in the variable x 
Then, the following holds 

1 



[We, 9) = 



W^a N (t,X)dtdX + Tu s{e N ). (41) 



£ (5d+2)/2 

Sketch of proof of Proposition 3.1. 

Using the short-hand notation Xs{t) '■= e~ QE *(l — x) {t/6\ we have 

{w e , 4>) =i/e [ 1 xs(t) (xs(H £ )S £ , U e (-t) 9e ) dt. (42) 
Je 

To compute the term U e (—t)(j) £ accurately, we use a projection over the overcomplete basis 
of L 2 (M. d ) obtained by using the so-called gaussian wave-packets, namely the family of functions 
indexed by (q,p) € K 2d defined by 



vIp&'O : = j^jdji ex P (~p- { x ~l)) ex P (- 

The point indeed is that, as proved by Combescure and Robert in [7], we have 
U £ (-t) cpl v (x,0 = O TuS (e N )+ 

(*-!)' 





i 


i exp ( 


-Pt- 

e 






exp^ 


if 

Jo 



, k-«tl 5 

TT I I 6Xp 



2e 

vl . 2, \ , qt-Pt-q-p 

— + n (q s ) jds 



in L°°([0, Ti]; L 2 (R d )). In other words, we have a quite explicit complex-phase/amplitude repre- 
sentation of the Schrodinger propagator when acting on the gaussian wave packets. 
This observation leads to writing, successively, in (42) 

( XS (H e )S e ,U e (-t) 9e ) = — / (x S (H e )S e ,<pl p )(<pl p ,U e (-t)<f> e )dqdp, 



(2ire) d 
1 



(xs(H £ )S £ ,ip qp ) (U £ (t)(p e q p ,(j) £ ) dqdp. 



Now, the idea is to replace the factor U £ (t)tp qp by its approximation derived above. Yet a few 
preliminary steps are in order. The first one uses the truncation in energy xs(H £ ), together 
with the functional calculus for pseudo-differential operators of Helffer and Robert (see [11] ), to 
replace this truncation by an explicit truncation near the set p 2 /2 + n 2 (q) — 0, up to small error 
terms. The second step consists in using the Parse val formula to write (we want to exploit the 
source term S £ on the Fourier side) 

for some function x(x) that truncates close to x — 0, and similarly 



(U e (t)<f$ tP ,<f> e ) = j—^ J x(y)j' m (U s (t)vl P ) (y)dydr), 
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These two steps explain the truncation factors \q and \i m the Proposition, which act close 
to the zero energy-level in phase-space (this is where functional calculus is used) and close to 
the origin in physical space. The last step consists in exploiting formula (43) in the obtained 
representation. 

Eventually, one obtains the desired formula. I 
3.3 Preparing for a stationary phase argument 

This slightly technical paragraph is devoted to proving that the obtained phase ijj in Proposi- 
tion 3.1 satisfies the assumptions of the stationary phase Theorem. 

Our main result in this paragraph is the Proposition after the following Lemma. 

Lemma 3.2. Let n 2 be any smooth refraction index. Then, the following holds 

(i) The stationary set associated with the phase ip in (37), defined as 

M x := {(t,X) = (t,q,p,x,£,y,ri) E [6,T X ] x R 6d s.t. V t ,x^(t, X) = and Imip{t, X) = 0} 
satisfies 

M x ={(t,q,p,x,£,y,r)) s.t. x = y = q = 0, £ = p, (t,p,r)) G M}, (44) 
where we recall that M = {(t,p,rj), X(t,0,p) — 0, £(t, 0,p) = 77, rj 2 /2 = n 2 (0)} by definition. 

(ii) We have, whenever m — (t,X) £ Mx, the relation 

Kei(D 2 ip lm ) = {{T,Q,P,X,Z,Y,H) e]0,+oo[xR 6d , X = Y = Q = 0, (45) 
E = P, rj T H = 0,B t (0,p)P + Tr} = 0, -H + D t (0,p) P + TVn 2 (0) = 0} . 

Note that this Lemma does not use the particular structure of our index. 

Proof of Lemma 3.2. A mere computation of Imi/j and V ip allows to write (44). Differenti- 
ating V ip once allows to write (45). For more details, the reader may check [-]. ■ 

With this Lemma at hand, our key result in this section is the following 

Proposition 3.3. Let n 2 be the refraction index defined in (21). We recall that the refocusing 
set M is computed in Lemma 2.4 and satisfies 

M = {(T R ,£,rj) s.t. £ = = (r, 0i, . . . , d -i), r = V2n 2 (0), |0i|<0 o }- 
Now, take any 

meMx = Ut,q,p,x,£,y,r)) s.t. x = y = q = 0,£,=p, 

(t,p, rj) E M, withp = (r, #i, . . . , Qd-i), and |0 X | < o j 

Then, we have 

Ker D 2 ij\ m =T m M x , 
where T m Mx denotes the space tangent to M x at point m. 

The remainder part of this subsection is devoted to the proof of Proposition 3.3. We begin 
by proving the Proposition in the case 

m = m := (T R ,0,p o ,0,p o ,0, -po), where p := (\/2n 2 (0), 0, . . . , 0) . 

We next generalize the result to other values of m, using the symmetries of the problem. 
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3.3.1 Proof of Proposition 3.3 when m = tuq 

The computation of T mo Mx on the one hand is rather easy 
Lemma 3.4. The space T mo Mx is given by 

T mo M x ={(T,Q,P,X,E,Y,H) s.t. X = Y = Q = T = 0, S = P = —H, P. Po = Q}. 

Proof of Lemma 3.4. This is a mere computation starting from the definition of the refocusing 
set M, as M = {(t,p,rj), X(t,0,p) = 0, E(t,0,p) = r), r/ 2 /2 = n 2 (0)}. ■ 

In order to determine Ker D 2 ifj\ mg the first step it to compute the matrices B t and D t involved 
in the linearized flow, see (3.3.2). 

Lemma 3.5. Let n 2 be the potential defined in (21). Then, we have 

D(T R ,0,p Q ) := ^(T R ,0,p ) = -I d , B(T R , 0,p ) := ^(T fl , 0,p Q ) = ^ ^ (46) 

where Id is the identity matrix, bu € M and Od-i is a square matrix of dimension d — 1 equal 
to 0. 

Proof of Lemma 3.5. 

We consider Xo(t, 0,p) = (xj^, 0,p), . . . ,Xq(£, 0,p)) the solution to (13) with initial data 
Xq(0, 0,p) = and x' (0, 0,p) — p. 

We recall that the index n 2 is invariant under the action of ©d,i(IR d ). Thus we first compute 
the components of D and B that are invariant under O^.i (R d ), namely their first column. We next 
compute the other columns by using the symmetries again, in conjunction with a perturbation 
argument. 

• Computation of §^(T R ,0,p o ) and §f (Tr,0,po) 
We start with ^(T R ,0,p ) for j > 2. We have 

ffs 4 s j (r Jl) o,( v ^(o) + E,o...,o)) -s J (r fl ,o,( v ^?(o),o...,o)) 

— HT R ,0,p Q ) = lim— ^ ^ ^ >-. 

Ot,i e^O £ 

Since the trajectory is radial we have 

Ej(T R ,0, (V2n 2 (0)+ £) 0...,0)) = H j (t R) 0, ( V^ 2 (0), . . . , 0)) -0, V j > 2. 

Hence, ^-(Tr,0,Po) = 0, V j > 2. A similar argument provides ^^-(T R ,0,p o ) = 0, V j > 2. 
There remains to determine the first coefficient of D. namely ^^(T R , 0,po). Since the trajectory 
is radial, and by conservation of the energy, we have for e small enough 

Si (T R ,0, (V2n 2 (0)+£,0,...,0)) =-(V2n 2 (0)+ £ ) , 

Si (r fl , 0, (V2n 2 (0),0,...,0)) = -^(0). 



Thus, 



(in := lim 

£->0+ 



5 (Tr, 0, (^2n2(0) + £, 0, . . . , 0)) - H (t r , 0, (^2^(0), 0, . . . , 0)) 
im — ^ '- ^ '- = -1. 
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• Computation of 85( ^°' Po) and ax( ^°' Po) (j > 2) 

Considering the symmetries of the problem, it is enough to consider the case j = 2: the other 
components may be determined using the same argument. 

We perturb the initial speed along the direction e%, by a factor e (see Figure 5). 




Figure 5: Perturbation of the initial speed 



Let X E (t) be the solution of the perturbed problem 

{x'J(t) = Vn 2 {X £ (t)), X e (0)=0, X' e (0)=p +ee 2 . 

We expand X e (t) with respect to e and obtain X £ (t) = X Q (t) + eXi(t) + . . .. With this notation 
we have X\{t) = §^(t) and X[(t) = §^(t)- To obtain the expansion in e, we go back to the 
previous case (J = 1) using a change of variables. Indeed, for e small enough, the trajectory is 
radial along the direction X' £ (0). Let (el, . . . , id) be a new basis defined by ij := O e Bj, with 



/cos(6> £ ) -sin(0 e ) 
sin(0 e ) cos(6» e ) 




Ia- 



d-2 







/ 



cos(6» e ) 



Pa 



sin 



Pl+e 2 ' 



\ 

Let X e be the coordinates of X s in (el, . . . , id). Since 0~ 1 Vn 2 (X e ) = Vn 2 (X e ), we clearly have 



X £ (t) = Vn 2 (X £ (t)), X e (0) = 0, X e (0) = ty^+pg.Q,..., 0)=p + O(£ 2 ). 
Hence it is clear that X £ (i) = Xo(t) + 0(e 2 ). Therefore, we recover 



X (t) + eXi(t) = O e (*o(t) + 0(e 2 ) = (J d + + O(e 2 ))(X (t) + 0(e 2 )), 



with 





f° 


Po 


... 0\ 




j_ 

Po 







... 


E := 









U 
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In other words, we have 

VteE, X (t)=X (t) and Xi(t) = EX (t). 
Since the Hamiltonian trajectory goes back to the origin at time T Rl we deduce 

QV 

— (T R ,0,p ) = Xx(T R ) = EX o (T R ,0, Po ) = ExO = 0. 
In the same way, we have 

§(T R APo) - X[(T R ) = EXq' (T R ) = ' IhPtirf px), . . .,x'q R (T r ) \ , 

o& \ Pa Pa J 

= t (0,l,0,...,0). 

The columns of B and D (for j > 3) are determined in the similar way. This leads to (46). ■ 

At this stage, we deduce the 
Corollary 3.6. KerD 2 tp\ m =T mo M X - 

Proof of Corollary 3.6. 

According to (45), we have 

Ker(D 2 ^ lm ) = {(T, Q, P, X, H, Y, H), X = Y = Q = 0, 

3 = P, r] T H = 0, B Tr (0, p) P + Tr) = 0, -H + D Tr (0, P )P + TVn 2 (0) = 0} . 

Since r\ = —po, we recover H = (0,i?2, • ■ ■ , Hd) ( m Cartesian coordinates). Since Vn 2 (0) = 0, 
we deduce that Dt r (0,p) P = H. According to Lemma 3.5, we deduce that H = —P. Finally, 
B Tr (o, p ) P = hence T = 0. Thus, 

KerD 2 i/> Uo = {(T, Q, P, X, E, Y, H), X = Y = Q = T = 0, P = S = -H, P.p = 0}. 

Using Lemma 3.4, the proof is complete. ■ 

3.3.2 Proof of Proposition 3.3 for any m 

In this subsection, we prove the 

o 

Lemma 3.7. Mm G M x , we have T m M x = KerD 2 i/>| m . 
Proof of Lemma 3.7. 

o 

The idea is to use a family of transformations which leave M x and n 2 invariant (in a sense 

o 

we define later), next to transport the equality Ker D 2 ip\ mQ = T mo Mx to any m G Mj. 

o 

Family of transformations. Let m — (t,q,p,x,£,y,rf) G Mx- We write m = 
(T Rl 0,p, 0,p, 0, -p) for some p G y/2n 2 (0) S^ 1 Thus, there exists R p G <0(R d ) such that R p (p) = 
p . We define the map R m : R 6d+1 — ► R 6d+1 by 

R m (t, q,p, x, £, y, 7?) = {t, Rp(q),Rp(y),R p (x),Rp(^),Rp(y),Rp(ri)) . (47) 
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By construction we have R m (to) = mo- 
Action on the tangent place. We have identified that the set Mx satisfies 
M x = {(t,q,p,x,£,y,rj), s.t. t = T R , q = x = y = 0, p = £ = -rj, p 2 /2 = n 2 (0)} 

s v— ' 

■-Mx 

n{ P = (r,e 1 ,...,e d _ 1 ) with |^| <e }. 

The set Mx is clearly invariant under the action of R m . Therefore, by restricting the domain in 

o o 

the variable 9\, it is clear that whenever m £ Mx, there exists a neighbourhood U of to in M x 

~~ o _ 

such that Uq := R m U C M x - Since the application R m is a linear map from U to Uq which 
satisfies R m (to) = too, we deduce 

Rrn (T m M X ) = T mo M X . 

Action on the kernel. We now compute the set i? m (Ker(Z) 2 -!/>| m )), as follows 

R\ n (Kev(D% m )) = {(T,R p Q,R p P,R p X,R p E,R p Y,R p H), s.t. X = Y = Q = 0, 
p.H = 0, B Tr (0, p )P + T p = 0, D TR (0 : p)P = H}, 

= {(T,Q,P,X,E,Y,H), s.t. X = Y = Q = 0, 

p.R~ 1 H = Q, B Ta (0 lP )R; 1 P + T P = 0, D TR {Q,p)R p 1 P = R- 1 H}. 

= {{T,Q,P,X,Z,Y,H), s.t. X = F = Q = 0, 

P o-H = 0, r p b Tr {o,p)r- 1 p + t Po = o, r p d Tr (o iP )r- 1 p = h}. 

On the other hand, we claim that 

R p B Tr (0, P )R- 1 = B Tr (0 iP0 ), R p D TR (0,p)R p - 1 = B Tr (0, Po ). (48) 
Assuming the above identity is proved, we immediately deduce 

R m (Kev(D%J)=KeiD 2 i; lmo . 

We conclude by writing 

Rm (Ker(D 2 tp\ m )) = KeiD% mo = T mo M x = R m (T m M x ). 

Thus, there only remains to prove (48). By construction of the potential we clearly have 

RpX(t,0,p) = X(t,0,p o ), as well as n 2 (R p x)) = n 2 (x) , 

whenever x/\x\ lies in the angular sector \6\\ < 9$. This provides 

D 2 n 2 D 2 n 2 D 2 n 2 

— (X(t,0,p )) = -^(R P X(t,0, P j) = R p —{X{t,Q,p))R-\ 

Therefore, using the differential equation (3.3.2) relating the time evolution of B t and D t , we 
recover the following system 

d 

l R p B{t,0,p)R p 1 = R p D(t,0,p)Rp 1 , R p B(0,0,p)R p 1 = Id, 



dt 

^ t R p D(t 1 0,p)R p 1 



D 2 n 2 

D 2 n 2 
Dx 2 



R 



D 2 n 2 
p ^x* 



(X(t,0,p))B(t,0,p)R- 



(R p X(t,Q,p))R p B(t,0,p)R p 1 
{X(t,0,p o ))R p B{t,0,p)R- 1 



R p D(0,0,p)R p 1 = 0. 
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Uniqueness of solutions to a differential system then gives 

Vfc, R p B t (0, P )R- 1 = B t (0,p ), R p D t (0,p)R; 1 - D t (0,p ). 
Relation (48) is proved. _ 

3.3.3 A useful byproduct of the proof of Proposition 3.3 

Lemma 3.8. Let n 2 be the refraction index defined in (21). Take any m £ Mx, written as 
m = (Tr, 0,p, 0,p, 0, — p) with p = ( -\J 2n 2 (0), 9\ , 62, • • • , Qd- 1 J according to Lemma 2.4- Then, 

tj}{rn) is constant on the set \9\\ < #o- 

Proof of Lemma 3.8. Considering the actual value of ift(m), various terms need to be consid- 
ered. The term / (p 2 /2 + n 2 (q e )) ds is clearly constant whenever \9i\ < 9 . The same statement 
Jo 

holds for the factor p t ■ q t . The only non-obvious factor is T t qt • Qt- As in the preceding proof we 
write 

T t (0,p)q t (0,p) -qt(P,p) =r t (0,p)q t (0,R; 1 p o ) ■ q t (0, R^po) 
= R p T t {0,p)R- 1 q t {0,p o ) ■ q t (0,p o ). 

There remains to write 

R p T t (0,p)R- 1 = R p (C t (0,p) + »D t (0,p)) • (A t (0, P ) + iB^p))- 1 R- 1 

= (R p C t (0,p)R- 1 +iR p D t (0,p)R p ~ 1 ) ■ (R p A t (0, P )R- 1 +iR p B t (0,p)R- i y 1 

= r t (o, Po ) 

for we already know that R p B t (0,p)R~ 1 = B t (0,po), R p D t (0,p)R p 1 = D t (0,po), and a similar 
proof establishes R p A t (0,p)R- 1 = B t (0,p ), R p C t (0,p)R p 1 = D t (0,p o ). u 

3.4 The stationary phase argument: Proof of item (hi) of our main 
Theorem 

The main result of the present section is 

Proposition 3.9. Let n 2 be the potential constructed according to 21. Select a source S G <S(R d ). 
Then, the following holds. 

(i) 7/supp f<S(£)^ n dle a = 0, we have 

V0eS(IR d ), (vT c -L e ,4>) = Tl)(5 (Vi), 

where (L £l (f>) is defined in (24) above (see also the Remark after Theorem 1.5), and dlg = {£ = 
(|£|, &i, . . . , 6d-i) such that Q\ = ±#0} (see definition 1.4). 

(ii) In the general case we have 

V0eS(]R d ), (uT e -L e ,4>) =o TuS (e°). 
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Proof of Proposition 3.9. Due to the fact that the stationary set Mx in the to-be-developped 
stationary phase argument has a boundary at 9\ = ±#o, the argument is in two steps. This is 
the reason why the above Proposition distinguishes between two cases. 

•• Proof of Proposition 3.9-part (i) 

Outside the stationary set Mx associated with the complex phase tp, the oscillatory inte- 
gral (41) defining (u^,</>) is of order 0(e°°). On the stationary set Mx and near the sup- 
port of a^v, the stationary set Mx is a submanifold without boundary, having codimension 
k = 6d + 1 — (d — 1) = 5d + 2. Indeed, thanks to the hypothesis on the support of S, we 
have suppajv H dMx = 0- 

Let us now come to the explicit application of the stationary phase Theorem to the oscil- 
latory integral (41). Writing p — (r, 9\, . . . , 9d-i) m hyperspherical coordinates, we define the 
application: 

7 : R Cld+1 n supp a N — ► M 5d+2 x S^ 1 

(t,q,p,x,(,y,rj) i — > (t - T R , q, x, y, £ - p, i] + p, r - ^2n 2 (0), 9 1; . . . , 9 d -i) 



The map 7 is a C°°-diffeomorphism between supp and 7 (supp ajv). Furthermore, we have 
by construction 

(i, X) £ M x n supp a N <^=> a = 0. 

The new coordinates (a, 9) are adapted to the stationary set Mx associated with ip. Making the 
change of variables (t,X) — 7 ~ 1 (a, 9) in the integral defining (uf e ,<j)) we have 

(w E ,<j>)=O s , Tl (e N )+ (49) 
1 



(5d+2)/2 



e 



~f(supp a N ) 



e^- 1 ^) (s{.W(.)P N ., ., -jj^J o 7 - 1 (o, d)x 3 (a, 9) r^ 1 dada(6), 



where da (9) denotes the standard euclidean surface measure on the unit sphere S d_1 , and X3 
is a truncation function on some compact set, a neighbourhood of Mx, whose precise value is 
irrelevant. Here we have used the non- stationary phase Theorem to reduce the original integral 
to an integral on a given compact set. 

Since for all point m £ Mx H suppajy we have Ker(_D 2 -0| m ) = T m Mx (Lemma 3.7), the 
function D 2 i\) is non-degenerate in the normal direction to Mx, which gives 



det( (0,0)) ^0. 



Furthermore, the projection of 7(supp a at) onto the space variable 9 is the angular sector 

IWe := {{9\, ... , 9d-i), 61 e] — So, #o[} , 

where lig denotes the projection (r, 6\, . . . , 9d-i) i-> (9%, . . . , We can now apply the station- 

ary phase Theorem in (49). Remembering that the codimension of the stationary set Mx associ- 
ated with ip is 5d + 2, we obtain that for any integer L there exists a sequence (Q2K^)W{o,....l} 
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of operators of order 21 such that 

(v£,4>) = 

exp (if sgn^^f^O, 
a ehQ de t^zi ( 0) e) 



exp -^o 7 - i (0,6») 



Q (.)S(.)4>*(.)Pn o 7-^3 (0, 0) d<r(0) 



exp(-^o 7 - 1 (O,0))^ £ £ Q 2 ,(9) ( (.,.,., 



t.=i 



.L+l 



O I £~ 1 ' SUp 

if<2L+d+3 



7 X3 



°7 _1 X3 J (O,0)d<7(0) 
+ 5>Tl ( £ ^) (50) 



: =I e + /7 £ + J77 e + 5)Tl (£ iV ), 
with the value 

d = (2vr)( 5rf + 2 )/ 2 (2n 2 (0))( d - 1 )/ 2 . 

The last line in ( ) serves as a definition of the three terms I e , II £ and IIL. and the L 00 -norm 
in 777 E is evaluated on a compact set of values of (a, 9), whose precise value is irrelevant. 

We compute these three contributions. Note that the retained value of the integer L remains 
to be determined at this stage. 

• Contribution of the remainder term 777 e in (50). 

This term is best studied by coming back to the original variables (t,X) instead of (a, 6). 
Expanding the fc-th order derivatives involved in this term, we clearly have 



777, = O ( e L+1 sup 

K<2L+d+3 



f x) [S{.)ct>*{.)P N (.,.,., -j. 



= O I e L+1 sup 

K<2L+d+3 



K 



(t,X) r N !•)•)• 



Hence, since 



P N (t, q,p, x) = ir-^ 4 det(A(t, q,p) + iB(t, q,p))~ 1/2 Q N (t, q,p, x), 



we recover 



777 e = O | e L+1 sup 

K<2L+d+3 



d (t, q ,p,y) {QN(t,q,p, (y-q t )/Ve)) 



Lastly, using (40) we have 



Efc 
^^~ J Pk,j(t,q,p,x), 

(k,j)ei N 



where pkj has at most degree k in x. We deduce 



777 E = Yl 



e t-J+ i + 1 sup 

K<2L+d+3 



d (t,q,p,y) {Pk,j{t,q,P, (y - qt)/Ve)) 



= E O {e*-S+ L+1 -*) = O (e L+1 -^ N -^) , 

(k,j)£l N 
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where we have used that j < 2N — 1 whenever (k,j) £ In (see (40)). There remains to chose 

L = 2N — 1 

to recover 

III £ = 0(e). 

• Contribution of II e in (50). 

This estimate is more delicate. Firstly, we have 



II, =Xy° 



Q 2e (d) S(.)^(.)P N o 7 -i X3 (0, 9) 



Hence, going back to the (t, X) variables again, and remembering that the relation (a, 9) = (0, 
implies y = q t = and t = Tr, we recover the identity 



II E = e * ° ( SU P 



K<2i 



U (t,q,p,y) 



y=q t =0,t=T R 



P N t,q,p, 



y-Qt 



where the L°°-norm is evaluated on some compact set of values of p. Now, inserting the exact 
value of P/v, we may write 



1=1 

L 



sup 

K<21 



U (t,q,P,y) 



E E 

i=i (kj)ei N 



(k,j)ei N 

e e s^~ j O [ sup 



y=q t =0,t=T R 



e 2 3 Pk,j t,q,p, 



K<21 



U (t,q,P,v) 



y=q t =0,t=T R 



Pk,j t,q,p, 



y -qt 
y - Qt 



Hence, using the fact that each pkj is a polynomial in its last argument, so that the above 
derivatives evaluated at y — q t — only leave the zero-th order term in the derived polynomial, 
we recover 



K £=l (k,j)£l 
L 



e e 2 J sup £ 

K<21 



-K/2 



o E E )=° E E =° £l/2 



<?=i (k,j)ei N 



1=1 (fc,i)G/ A 



where we have used that k — 2j > 1 whenever (fc, j) G J 



• Contribution of 7 £ in (50). 

The integral defining I £ has the following more explicit value, where p = (^/2n 2 (0), 9\ , . . . , 0<j_i) 
namely 



exp(( -r/>(TR,0,;p,0,p,0,-p) 



n e /, det(^F i (0,e)) 

det(A(T R ,0,p) + iS^.Cp))- 1 / 2 5(p) £*(-p) d9 x . . . d9 d . u 
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On top of that, we have 

i/j(T R ,0,P,0,P,0-p) = £ R ( IPs( ° 2 P)|2 +n 2 (g,(0,p))) ds t 

while the fact that n 2 is radial implies that ip(T R , 0,p, 0,p, 0, — p) = ip{T R , 0,Po, 0,poj 0, — Po) 
whenever p E Ig . For the same reason, we also have whenever 9 € Hele the relation 



det(^^(O,0)) det (^=1(0,0)) 
together with the identity, valid whenever p G Ig Q , 

det(A(T R ,0,p) + iB(T R ,0,p))-^ 2 ^det(A(T R ,0, Po ) + iB(T R ,0,p Q ))- 1 ^. 
Eventually, we have obtained 

(i [ Tr f\ P s(o, Po y 2 



. , ., +« (9.(0,Pto)) J ds\ 

I £ = C n2 , d e\ £ Jo V M %)0*(-p)^ o (p), (51) 

with 

C TR .<i == — T^r, ; v det(A(T R ,0,p Q ) + iB(T R ,0,p o ))-^ 2 . 

det (^=i(o,0)J 

This ends the proof of Proposition 3.9-part (i). 
•• Proof of Proposition 3.9-part (ii) 

In that case, the argument is essentially the same (a stationary phase argument in the vari- 
able a) , up to a convenient use of the dominated convergence Theorem (to deal with the variable 
01, and more specifically with the boundary 9\ — ±9q). 

Namely, we first write, as in the proof of part (i) of the Proposition, 

(w e ,<f>)=Os,T x (z N )+ (52) 



(5d+2)/2 



e 



e 

j(supp a N ) 



i+^M (s(.)H)P N (.,.,. ,-j=Y) o 1 - 1 (a,9) X 3(a,9)r d - 1 dada(9), 



where X3 is a truncation function on some compact set, a neighbourhood of Mx, whose precise 
value is irrelevant. Here we have used the non-stationary phase Theorem to reduce the original 
integral to an integral on a given compact set. The key point now lies in writing, 

(%,</>) =O s ,Ti(e rr )+ (53) 
d ^ 6 ) J daei^ 1 ^ (s(.W(.)P N (•,-,.,-^)) o 7 -\a,6) X3 (a,0)r d -^ . 

v v ' 

With this formulation in mind, our next objective is to prove that whenever 77 > is a small 
parameter we have 

f da(9)\J E (9)\<C V7 (54) 

J\0i±e o \<v 
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for some C > independent of e and rj. It is clear indeed that the upper-bound (54), in 
conjunction with part (i) of the Proposition, provides a complete proof of Proposition 3.9-part (ii). 

Let us now concentraate on the case \0\ — 0o| < rj (the proof in the case |0i + 0o| < 77 is the 
same). 

In order to prove (54), we fix a value (0°, . . . , ano - we P rove that, given (0°, . . . , 02_i), 

there is an 77 > 0, and a C > independent of e, such that 

V0 such that |0 - (0 O .0°, . . . , < »?, we have |J e (0)| < C. (55) 

Covering the whole set {0 € § d_1 ; |0i — 9q\ < ??} by finitely many sets of the form {|0 — 
(00-02! • • • j ) I < ^} clearly provides the desired relation (54) once (55) is proved. 

Now, relation (55) results from an application of the stationary phase Theorem, with complex 
phase and with parameter. Here a is the variable used for the stationary phase itself, while 
is the parameter, and ip o 7 -1 is the complex phase. We introduce the short-hand notation 
0° = (00: (0')°) = (@0i 0°) • • • ) @d-i) f° r convenience. It has already been established 6 that 

ImfV-c^- 1 ) (a, 6) > 0, V(a,0), 
Im (V) o 7- 1 ) (a = 0, = 0°) = 0, 
V Q (V>°7 -1 ) (a = O,0 = 0°) = O, 

det( D ^)(a = 0,0 = 0^0 

Therefore, the stationary phase theorem with parameter ensures that close to = 0° there is an 
expansion of the form 



j e (0) = e *W/« j^Q 2 ,(a Q ) ^(.)^(.)P N o 7 - 1 X3 (.)jj (0)j+i?(e,i 



for some smooth functions <f> and R(e, L, 0), where the Q21S are differential operators of order 11 
in the variable a, and, for any function u(a, 0), the notation u°(0) refers to any smooth function 
m°(0) that belongs to the same residue class than the original function u(a,9) modulo the ideal 
generated by V q V ;o 7 _1 (q!, 0) (see Hormander [12], sect. 7.7, for the details). With this notation, 
we actually have <p = (1/1 o Besides, the remainder term R satisfies as the term III £ in 

the previous step an estimate of the form 



\R{e,L,0)\ <C L s L+1 sup 

\K<2(L+1) 



^((§(.)?(-)^(-,-,-, 77 |))°7- 1 X3(.) 



for some constant Cl > independent of e, and provided is close to 0° (independently of e). 
These two ingredients immediately provide, using the same estimates as we did for the terms 
III e and II e above, the upper-bound, valid for close to 0°, 



\J e (e)\<c (Xy ^g 2 ,(9 Q )^(.)^(.)p A r^.,.,-^Jo 7 - 1 X3 (.)jj (0)J+i?( £ ,L, 

Gathering powers of e as in the previous part of the proof, provides the upper bound 

\U6)\<C, 



6 Stricto sensu, these relations have only be proved when < 80, and we here extend the result to the case 
9i = do- This is allowed due to the invariance of the phase on the parameter 8 whenever \8±\ < 8q - Lemma 3.8. 
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where C does not depend on e and 9 is close to 9°, independently of e. Point (55) is proved. 

We immediately deduce that (54) holds, and the proof of Proposition 3.9 - part (ii) is com- 
plete. 



3.5 Conclusion 

Gathering the intermediate result in Proposition 1.6, together with Proposition 3.9, gives item 
(iii) of Theorem 1.5, by conveniently choosing the parameters 6, 9, Tq and T%. 
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